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1 Introduction 

The discovery of a standard model-like Higgs boson at the CERN LHC (Tj-Hl opens a new phase 
in the understanding of the standard model (SM) of particle physics. The search for additional 
Higgs-like particles and the measurement of their properties provide complementary ways to 
test the validity of the SM and to test for the presence of physics beyond it. 

This analysis describes a search for new resonances in the diphoton invariant mass spectrum, 
using data corresponding to an integrated luminosity of 19.7 fb 1 collected with the CMS de¬ 
tector at a center-of-mass energy of 8 TeV at the CERN LHC. Despite the large nonresonant 
background, the diphoton decay mode provides a clean final-state topology that allows the 
mass of the decaying object to be reconstructed with high precision, exploiting the excellent 
performance of the electromagnetic calorimeter of the CMS experiment. The analysis searches 
for local excesses that could be due to the production of particles that decay into two photons 
with mass in the range from 150 to 850 GeV. Both narrow and wide resonances are investi¬ 
gated with natural widths ranging from 100 MeV to 10% of the resonance mass. This search 
covers the diphoton mass range above that investigated in |2j[4l. The ATLAS experiment re¬ 
cently published a similar search for narrow resonances in the diphoton final state in the mass 
range between 65 and 600 GeV at a center-of-mass energy of 8 TeV H5j. Previous searches for 
resonant diphoton processes have been performed at the Tevatron by DO U 6 j and CDF [7j at a 
center-of-mass energy of 1.96 TeV and by the ATLAS (8;1 and CMS (U experiments at the LHC 
at a center-of-mass energy of 7 TeV. 

Several models of physics beyond the SM, such as the two-Higgs-doublet model (2HDM) IflOl . 
motivate the search for additional high-mass resonances in the diphoton channel. Generally, 
these models provide an extension of the Higgs sector, where a total of five Higgs bosons are 
predicted by the theory. The mass spectrum of the 2HDM can be split into two regions: a light 
SM-like Higgs boson h with mass around 125 GeV and the remaining physical Higgs bosons, H, 
a scalar. A, a pseudoscalar, and H*, clustered at an equal or higher scale with niu ~ » 6 \ ~ m H +- 
Under the assumption that the newly observed Higgs boson is the light CP-even Higgs scalar 
of the 2HDM, the consistency of its couplings with those predicted by the SM pushes the model 
close to the so called alignment limit Ifffl . where certain decay modes of heavy neutral Higgs 
bosons vanish, including H —> VV (where V is a vector boson), H —> hh, and A —> Zh. At the 
same time, decays of H and A to 77 and tt become increasingly important and the electroweak 
production modes, such as vector boson fusion or production in association with a W or a Z 
boson, are predicted to be suppressed. Therefore the production of both H and A is dominated 
by gluon fusion. The absence of tree-level flavor-changing neutral currents in multiple-Higgs- 
doublet theories is guaranteed by the Glashow-Weinberg condition This condition is 

satisfied in the 2HDM by four discrete combinations of the Yukawa couplings of the fermions 
to the Higgs doublets. In the Type I scenario all fermions couple to one doublet, while in Type II 
up-type quarks couple to one doublet and down-type quarks and leptons couple to the other. 
A detailed description of other scenarios is given in Ref. ItTOj . 

Given the general character of this search, the results can also be interpreted in terms of differ¬ 
ent spin hypotheses for the new particle. The Landau-Yang theorem flT3lfL4l forbids the direct 
decay of a spin-1 particle into a pair of photons. The cases of spin-0 and spin-2 are investigated 
in this analysis. Spin-2 particles decaying into two photons are predicted by other extensions 
of the SM such as the Randall-Sundrum Ifl5l and the Arkani-Hamed-Dimopoulos-Dvali Ifl 6 l 
models. These theories predict a distinct higher-dimensional scenario, which provides an ap¬ 
proach to the hierarchy problem alternative to supersymmetry. The particle predicted in this 
context, the graviton, can have a mass in the TeV range and thus be observed at the LHC pref- 
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2 The CMS detector 


erentially in its decay into two gauge bosons, such as photons. 

The paper is organized as follows: after a brief description of the CMS detector in Section [2j 
the data and simulated samples are described in Section |3j while the reconstruction and iden¬ 
tification of photons is detailed in Section [4] The diphoton vertex identification is covered 
in Sectio n |4.3[ followed by the description of the event selection and classification in Section [5] 
Sections [Oja nd |6.2| describe, respectively, the signal and background models used for the inter¬ 
pretation of the data, and Section 6.3 discusses the associated systematic uncertainties. Finally, 
in Section[7]the model independent results of the search for new diphoton resonances and their 
interpretation in terms of the standalone production and decay rates for H and A within the 
2HDM are discussed. We express these results in terms of the appropriate 2HDM parameters. 


2 The CMS detector 

The central feature of the CMS detector is a superconducting solenoid, of 6 m internal diame¬ 
ter, providing an axial magnetic field of 3.8 T along the beam direction. Within the field volume 
there are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter 
(ECAL), and a brass and scintillator hadron calorimeter (HCAL). Charged particle trajectories 
are measured by the tracker system, covering 0 < (p < 2n in azimuth and \rj\ < 2.5 in pseudo¬ 
rapidity. 

Muons are measured in gas-ionization detectors embedded in the steel return yoke. The ECAL, 
which surrounds the tracker volume, consists of 75848 lead tungstate crystals that provide 
coverage in pseudorapidity \rj\ < 1.48 in the barrel region (EB) and 1.48 < \rj\ < 3.00 in two 
endcap regions (EE). The EB modules are arranged in projective towers. A preshower detector 
consisting of two planes of silicon sensors interleaved with a total of 3 radiation lengths of lead 
is located in front of the EE. In the region \t]\ < 1.74, the HCAL cells have widths of 0.087 in 
pseudorapidity and azimuth. In the plane, and for )/| < 1.48, the HCAL cells map onto 
5x5 ECAL crystal arrays to form calorimeter towers projecting radially outward from close to 
the nominal interaction point. At larger values of \rj\, the size of the towers increases and the 
matching ECAL arrays contain fewer crystals. Within each tower, the energy deposits in ECAL 
and HCAL cells are summed to define the calorimeter tower energies, subsequently used to 
calculate the energies and directions of hadronic jets. A more detailed description of the CMS 
detector, together with a definition of the coordinate system used and the relevant kinematic 
variables, can be found in Ref. lfT7l . 

Photons generally deposit their energy in a group of crystals of the ECAL called "cluster". Re¬ 
construction of the photons used in this analysis is described in Section|4j and uses a clustering 
of the energy recorded in the ECAL, known as a "supercluster", which may be extended in 
the (p direction to form an extended cluster or group of clusters Ifl 8 | . Several procedures are 
used to calibrate the energy response of individual crystals before the clustering steps The 
changes in the transparency of the ECAL crystals due to irradiation during the LHC running 
periods and their subsequent recovery are monitored continuously, and corrected for, using 
light injected from a laser system. The calibration of the ECAL is achieved by exploiting the (p 
symmetry of the energy flow, and by using photons from rr° —> 77 and ?/ —> 77 decays and 
electrons from W — > ev e and Z — > e ! e decays f!9| . 
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3 Data samples and simulated events 


The events used in the analysis are selected by two diphoton triggers with asymmetric trans¬ 
verse momentum thresholds (pj), 26 and 18GeV or 36 and 22GeV on the leading and sub¬ 
leading photons respectively, depending on the data taking period, and complementary pho¬ 
ton selections. One selection requires a loose calorimetric identification based on the shape of 
the electromagnetic shower and loose isolation requirements on the photon candidates, while 
the other requires only that the photon candidates have a high value of the Rg shower shape 
variable. The R 9 variable is defined as the energy sum of 3 x 3 crystals centered on the most 
energetic crystal in the supercluster divided by the energy of the supercluster. Photons that 
convert before reaching the calorimeter tend to have wider showers and lower values of R 9 
than unconverted photons. High trigger efficiency is maintained by allowing both photons to 
satisfy either selection. The measured trigger efficiency is above 99.8% for events satisfying the 
diphoton preselection discussed in Section 4.1 required for events entering the analysis. 


Monte Carlo (MC) signal and background events are generated using a combination of pro¬ 
grams. Full detector response is simulated with GEANT4 [(20| . Multiple simultaneous pp in¬ 
teractions (pileup) are simulated, and the events are weighted to reproduce the distribution of 
primary vertices observed in data. The interactions used to simulate pileup are generated with 
the same version of PYTHIA l2ll . v6.424, that is used for other purposes as described below. 
The PYTHIA tune used for the underlying event activity is Z2* If22l . New resonances X are sim¬ 
ulated with a natural width of 0.1 GeV which is smaller than the value of the mass resolution 
in the energy range considered. Both spin-0 and spin-2 hypotheses are considered. Simulated 
signal samples of X —> 77 events are generated with PYTHIA ETl for the gluon fusion process 
with the following mass hypotheses: 150, 200, 250, 300, 400, 600 and 850 GeV. Interference 
between the signal and the background is studied using SHERPA 2.1.0 Il23li24 j. The exclusion 
limits detailed in Section [7] are computed using either a signal plus interference template given 
by SHERPA or the Breit-Wigner theoretical model described in Section [6.11 The inclusion of the 
interference in our signal model affects the expected and observed upper limits on the signal 
yields at the level of 3% or less. The effect is therefore negligible compared to the systematic 
uncertainties, described in Section [63] and is not included in the final results. Simulated back¬ 
grounds include the diphoton continuum process involving two prompt photons, generated 
with SHERPA 1.4.2 H23H . and processes where one of the photon candidates arises from misiden- 
tified jet fragments, simulated with PYTHIA. These two contributions represent 97% of the total 
background in this analysis. A less than 3% contribution is expected from QCD events in which 
both photons candidates arise from misidentified jet fragments; this contribution is neglected 
in the analysis. 


4 Photon reconstruction and identification 

Photon candidates for the analysis are reconstructed from energy deposits in the ECAL using 
algorithms that constrain the superclusters in f/ and (p to the shapes expected from photons 
with high pj. The clustering algorithms account for about 98% of the energy of the photons, 
including those that undergo conversion and bremsstrahlung in the material in front of the 
ECAL. Groups of clusters are used to form superclusters. In the barrel region, superclusters are 
formed from five-crystal-wide strips in ?/, centered on the locally most energetic crystal (seed), 
and have a variable extension in cp to take into account the effect of the magnetic field on elec¬ 
trons from photons showering before ECAL. In the endcaps, where the crystals are arranged 
according to an x-y rather than an rj-cp geometry, matrices of 5 x 5 crystals, which may par¬ 
tially overlap and are centered on the seed crystal, are summed if they lie within a narrow cp 
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4 Photon reconstruction and identification 


road. About half of the photons convert into e ' e pairs in the material upstream of the ECAL. 
If the resulting charged-particle tracks originate sufficiently close to the interaction point to 
pass through three or more tracking layers, conversion track pairs may be reconstructed and 
matched to the photon candidate. 

The energy containment of the photon showers in the clustered crystals and the shower losses 
due to conversions in the material upstream of the calorimeter are corrected using a multivari¬ 
ate regression technique based on a boosted decision tree (BDT) 11251 1261 , which uses as input a 
collection of shower shape and kinematic variables, together with the energy measured in the 
preshower for events with photons in the endcaps. Corrections are derived from simulation. In 
order to correct for residual discrepancies between simulation and data, the MC simulation is 
tuned to match the energy resolution observed in data, while data are calibrated to match the 
energy response in simulation. More details about photon reconstruction can be found in Ifl8l . 

The photon candidates used in the analysis are required to be within the fiducial region, \t]\ < 
2.5, excluding the barrel-endcap transition region 1.44 < \rj\ < 1.57. 

4.1 Photon preselection 

The photon candidates in this analysis are required to satisfy preselection criteria similar to the 
trigger requirements: 

• Pt > 33 GeV and p^. 2 > 25 GeV, where Pj' and p'-f 2 are the transverse momenta of 
the leading and subleading photons, respectively; 

• a selection on the hadronic leakage of the shower, measured as the ratio of energy in 
HCAL cells behind the supercluster to the energy in the supercluster; 

• a photon identification based on isolation and shape of the shower with lower thresh¬ 
olds with respect the one used in the final photon selection described in Section |4~2j 

• an electron veto, which removes the photon candidate if its supercluster is matched 
to an electron track with no missing hits in the innermost tracker layers. 

The efficiency of the photon preselection, measured in data using a "tag-and-probe" tech¬ 
nique If27l . ranges from 94% to 99% 0. The efficiency of all preselection criteria, except the 
electron veto requirement, is measured using Z —> e 1 e events, while the efficiency for pho¬ 
tons to satisfy the electron veto requirement is measured using Z —> ppy events, in which the 
photon is produced by final state radiation, which provide a more than 99% pure source of 
prompt photons. 


4.2 Photon identification and selection 


Photon identification is performed by applying selection requirements on a set of discrimi¬ 
nating variables. In this analysis, the selection is optimized separately in four different pseu¬ 
dorapidity and R 9 regions 0. The variables used to suppress the background due to the 
misidentification of jets with high electromagnetic content are 


• the sum of the transverse mome nta o f all the tracks coming from the vertex chosen 
for the event (described in Section 4.3' within a veto cone of A R = \/ (A //) 2 + (A tp) 1 = 
0.3 around the photon candidate's direction; 


• the sum of the ECAL energy deposits in crystals located within a veto cone of AR = 
0.3 around the supercluster position, excluding the photon; 


• the sum of the energies of HCAL towers whose centers lie within an annular veto 
region of outer radius AR 0 = 0.4 and inner radius AR, = 0.15, centered on the ECAL 
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supercluster position; 

• the ratio between the sum of HCAL tower energies within a veto cone of size A R < 
0.15 centered on the ECAL supercluster position, and the energy of the supercluster; 

• the spread in // of the electromagnetic cluster, computed with logarithmic weights 
and defined as 


A = Ea, 'fy ~ g where 


VV 


L™i 


V = 


L^iVi 

L^i 


and Wi = max 0,4.7 + In 


£5 


x5 


( 1 ) 


and the sum runs over the 5x5 crystal matrix around the most energetic crystal 
in the supercluster. E 5 x 5 is the energy of the 5x5 crystal matrix and //,• = 0.0174//,, 
where //, is the ?/ index of the zth crystal. This variable represents the second moment 
of the energy distribution along the // coordinate. 

A more detailed description of the isolation requirements can be found in Ifl 8 l . The energy 
deposited within the isolation cones is corrected for the contribution from pileup and the un¬ 
derlying event using the FASTjET technique |28| . 


4.3 Diphoton vertex identification 

The diphoton mass resolution has contributions from the resolution of the measurement of the 
photon energies and the angle between the two photons. The opening angle resolution strongly 
depends on the determination of the interaction point where the two photons were produced. 
If the vertex from which the photons originate is known with a precision better than 10 mm, 
the experimental resolution on the angle between them makes a negligible contribution to the 
mass resolution. 

No charged particle tracks result from photons that do not convert, so the diphoton vertex is 
identified indirectly, using the kinematic properties of the diphoton system and its correlations 
with the kinematic properties of the recoiling tracks. If either of the photons converts, the 
direction of the resulting tracks can provide additional information f 2 ]|. 

The efficiency for finding the correct vertex for a diphoton resonance of mass above 150 GeV is 
between 79% and 92% and increases with the mass of the resonance. 


5 Event selection and classification 

The analysis uses events with two photon candidates satisfying the preselection and identifi¬ 
cation requirements, and with pj 1 > m 77 /3 and p'j 2 > m 77 / 4. The use of pj thresholds scaled 
by niyy prevents a distortion of the low end of the m 77 spectrum that would result from fixed 
thresholds l29l . This strategy has two other main advantages. Background modeling is sim¬ 
plified because the shapes of the invariant mass distributions in the different event classes are 
similar after this selection. In addition scaled pj thresholds allow for tighter selection criteria 
with the mass increase, while preserving the acceptance of the selection. 

The search sensitivity is increased by subdividing the events into classes, according to indica¬ 
tors of mass resolution and predicted signal-to-background ratio. Two simple classifiers are 
used: the minimum Rg and the maximum pseudorapidity of the two photons. Photons with a 
high value of the R 9 variable are predominantly unconverted, have a better energy resolution 
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5 Event selection and classification 


Table 1: Definition of diphoton event classes 

Class tj criterion Rg criterion 

0 max(|//|) < 1.44 min(R 9 ) > 0.94 

1 max(|//|) < 1.44 min(R 9 ) < 0.94 

2 1.57 < max(|^|) < 2.50 min(R 9 ) > 0.94 

3 1.57 < max(j^j) < 2.50 min(R 9 ) < 0.94 
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Figure 1: Signal efficiency times acceptance as a function of the mass for a spin-0 scalar res¬ 
onance, produced via gluon-gluon fusion, with natural width equal to 0.1 GeV. The shaded 
region indicates the systematic uncertainty. 


than those with a lower value, and are less likely to arise from misidentification of jet frag¬ 
ments. Similarly, photons in the barrel have better energy resolution than those in the endcap 
and are less likely to be incorrectly identified. The classification scheme groups together events 
with good diphoton mass resolution, resulting from photons with good energy resolution and 
with better signal-to-background ratio. The event class definitions are shown in Table [lj 


Figure [l] shows the signal efficiency times the acceptance of the event selection as a function 
of the mass hypothesis for a narrow spin-0 scalar resonance produced via gluon-gluon fusion. 
The shaded region shows the systematic uncertainty due to different sources as described in 
The acceptance over the full mass range is also evaluated for the pseudoscalar 


Section 6.3 


hypothesis and found to be compatible with that measured for the scalar resonances within 
the uncertainties. For the spin-2 scenario the corresponding numbers are about 5-10% smaller, 
because of differences in acceptance between the two models. Figure [2] shows the diphoton 
invariant mass distribution for selected events in data and MC simulation, normalized to an 
integrated luminosity of 19.7 fb -1 , for all event classes combined. The bin width of this distri¬ 
bution is chosen to be narrow enough to properly display wide resonances. However a data 
driven technique is exploited in this analysis for the estimation of the background, as detailed 
in Section [6721 
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Figure 2: Diphoton invariant mass distribution for the selected events in data and simulation. 
Background processes are represented by the filled histograms. The shaded band represents 
the Poisson uncertainty in the MC prediction.The prompt diphoton events and the photon 
plus jets events are shown separately The ratio between data and MC is displayed bin by bin 
in the bottom plot. In this figure, the leading-order cross sections for the background processes 
are scaled to next-to-leading-order (NLO) predictions with scale factors derived from CMS 
measurements at 7 TeV I30ll3l11 . 


6 Statistical methodology 


In order to assess the compatibility of the data with the presence of a diphoton resonance, we 
make a hypothesis test based on a frequentist construction 1321133 1. An unbinned maximum- 
likelihood fit in a sliding window range to the diphoton invariant mass distributions in all the 
event classes is made using a parametric model for the signal and a background shape obtained 
directly from data. A scan of the signal mass and the signal width is performed. The signal and 
background normalizations are allowed to float, together with the parameters which describe 
the background shape. The log-likelihood ratio is used as test statistic. The signal model is de¬ 
rived from MC simulation as described in Section [O] with corrections determined from data- 
MC comparisons applied. The functional form of the background distribution is determined 
by fitting the measured m 77 distribution as detailed in Section 6.2 Systematic uncertainties de¬ 


scribed in Section[63|are incorporated into the analysis via nuisance parameters and are treated 
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6 Statistical methodology 


according to the frequentist paradigm j34j|. 


6.1 Signal parametrization 


To construct the signal fitting function from simulated events, the reconstructed rn r - f distribu¬ 
tion of each event class is described in terms of a parametric function of the hypothetical signal 
mass nix- This procedure can be extended in a simple fashion to allow for an additional free 
parameter, the natural width Tx of the new resonance, by convolving a resolution function, 
which takes into account the detector response, with a theoretical line shape. The convolution 
of the detector response, described in this section, with the theoretical lineshape, described in 
section 6.1.1 is applied for any natural width hypothesis. 


The detector response is parametrized in terms of the relative difference between the recon¬ 
structed diphoton mass m reco and the true mass mtrue, } l = (»beco — «hrue) /»6rue- The resolution 
function R is obtained by fitting the response distribution with an analytic function, namely 
the sum of two single-sided Crystal Ball (CB) functions ll35l with common mean }Iq and width 
(7, and different values of n and oc. The Crystal Ball function combines a Gaussian core and a 
power-law tail with an exponent n to account for incomplete photon energy containment in the 
cluster related to the material in front of the calorimeter, and for other reconstruction effects: 


/cb(m) = 


( —tL 

) 

I N 

V V / 27T< 


(h)”«p(-¥)(h 


CL 


F~Fo \~ n 
cr ) 


} l ~Vo 

cr 


> or, 


ICTo < a . 


( 2 ) 


The parameter a defines the transition between the Gaussian and power-law functions. The fit 
to the response distribution in the first class of events for a simulated signal with nix = 150 GeV 
is shown in Fig. [3] The resolution in \i (i.e., the a in Eq. (|2]l) improves by roughly 20%, from 


8 TeV 




Figure 3: Fit to the detector response distribution for nix = 150 GeV with two single-sided 
Crystal Ball functions (solid curve), displayed with linear (left) and logarithmic (right) scales. 
The dashed and dotted curves show the individual Crystal Ball components. 


0.021 to 0.016, when the resonance mass increases from 150 to 850 GeV. The resolution function 
R is constructed so that u depends continuously on fn truc via a quadratic polynomial function. 
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6.1.1 Parametric line shape 

The theoretical line shape of unstable particles is modeled with a Breit-Wigner (BW) distribu¬ 
tion with a mass-dependent width l36l : 


g B wM m x,r x ) 


N 

(;m 2 — mj ,) 2 + 


(3) 


The effect of the proton parton distribution functions (PDF) on the signal shape for a wide 
resonance (with r x = 0.1M X ) is investigated with high-mass Higgs boson samples produced 
with the POWHEG generator Il37l . The shape obtained after convolution with the PDF is still 
well described by a BW function within an accuracy better than 1%. 

The final signal model is obtained from the convolution of the BW in Eq. ([3]) and the response 
function R: 

§s( m ry\ m x, T x ) = R(m 77 \m) ®gBwMm x ,r X ). (4) 

This signal model is validated by fitting the reconstructed mass distribution for a simulated 
signal with T x = 0.1 GeV, as shown in Fig. [4] 

8 TeV 8 TeV 




Figure 4: Parametrized signal shape for a signal with ?n x = 150 GeV and F x = 0.1 GeV, dis¬ 
played with linear (left) and logarithmic (right) scales. The solid curve shows the result of a fit 
to the simulated events (points with error bars); the dashed curve represents the Breit-Wigner 
component of the model. 


6.2 Background modeling 

The modeling of the background relies entirely on the data, so there are no systematic un¬ 
certainties due to potential mismodeling of the background processes by the MC simulation. 
Because the exact functional form of the background in each event class is not known, the para¬ 
metric model must be flexible enough to describe a variety of potential underlying functions. 
Using an incorrect background model can lead to biases in the measured signal yield, which 
can strongly reduce the sensitivity of the analysis to any potential signal. The procedure used to 
determine the background fitting function, which results in a negligible bias, is presented here. 
In this study the bias on the fitted signal yield is defined as the difference between the number 
of fitted signal events and the number of expected signal events. A set of analytical functions 
that could describe the unknown true background distribution in data is first determined. The 
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6 Statistical methodology 


five functions considered as possible truth models are analytical forms that are used in dijet 
resonance searches H38l to describe both data and QCD predictions and models that are fre¬ 
quently used in diphoton resonance searches EJ • The bias study procedure described below 
is performed considering each function among the available set of models. The candidate fit¬ 
ting function that performs best in this study is fo, the product of exponential and power-law 
functions. 

/o(m) = e“ Pim m“ P2 . (5) 

As an example, the fits of this model to the data in the four event classes are shown in Fig. [5] 
for the fit range [240, 640] GeV, used for searching for a peak near 350 GeV. In order to check 


19.7 fb' 1 (8 TeV) 



19.7 fb" 1 (8 TeV) 




19.7 ft)- 1 (8 TeV) 



Figure 5: Fits to the diphoton mass distributions in the four classes of events in the window 
[240, 640] GeV chosen for searching for a peak near 350 GeV using the fo model and assuming 
no signal. The ratio between data and fo is displayed bin by bin in the bottom plot. 


that a background fit model results in a negligible bias in the fitted signal yield, we construct 
pseudo-data sets with the other four analytical models by randomly drawing diphoton mass 
values from them. The number of background events generated in each set is equal to the 
number of events observed in data in a fixed mass range. Background-only pseudo-data sets 
are then fitted with the signal + background probability density function. The criterion for 
the bias to be negligible is that it must be five times smaller than the statistical uncertainty in 
the number of fitted signal events for all four types of pseudo-data sets across the entire mass 
region of interest. When this criterion is satisfied, any potential bias from the background fit 
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function can be neglected in comparison with the statistical uncertainty from the finite data 
sample. The functional form of the truth model being fixed, the lower and upper bounds of the 
fit range are varied for each hypothetical resonance mass in order to minimize the bias. The 
desired fit range for each hypothetical mass is the one that gives the smallest bias for all the 
truth models considered. This study is performed for a resonance with width equal to 10% 
of its mass, as the bias tends to increase with increasing resonance width. The optimal lower 
and upper bounds of the fit range are parametrized as a function of the resonance mass with 
polynomial functions. 

As a closure test, the obtained fit ranges are used to compute the bias on the number of fitted 
signal events as a function of the assumed resonance mass. We find that the fo function pro¬ 
duces a sufficiently small bias for all four truth models in all event classes and for any width of 
the resonance up to 10% of the mass. We therefore use this background function to fit the data. 

Our approach to extract the numbers of signal and background events cannot be used above 
= 850 GeV because of the very small number of events in the data. This is therefore the 
highest value of mass considered in this search for new resonances. The lowest value of m\ 
considered is 150 GeV, and the fitted range in m 77 is [130,1000] GeV. The maximum value of 
the width of the resonance is fixed to 10% of the resonance mass itself. This value is limited by 
the width of the resonance mass points at the edges of the range (150 and 850 GeV) which have 
to be included in the fitting range ([130,1000] GeV) within at least one sigma. 

6.3 Systematic uncertainties 

The experimental systematic uncertainties can be separated into those related to the yield and 
those related to the signal shape. A log-normal prior is assumed for the uncertainties in the 
class yields, while the shape uncertainties are incorporated as parametric variations of the 
model. 

The normalization uncertainty related to the integrated luminosity is 2.6% Il39l . The photon- 
related uncertainties are the same as in Q. The uncertainty in the energy scale is conservatively 
increased to 0.5% in the barrel and to 2% in the endcaps to take into account additional non- 
linearities. Systematic uncertainties related to individual photons are then propagated to the 
signal model, where they result in uncertainties in the peak position and width. The 1% (2.6%) 
normalization uncertainty in the barrel (endcap) related to the offline photon identification is 
taken from the largest uncertainty in the data/MC scale factors computed with Z —» e 1 e 
events using the tag-and-probe technique. A 1% normalization uncertainty is also assumed 
in the trigger efficiency A normalization uncertainty of 2-5% is included to take into account 
possible event class migration due to R 9 selection efficiency uncertainties in both barrel and 
endcap. More details about trigger and R 9 efficiencies can be found in Ref. 0. 

Systematic uncertainties estimated for the SM-like Higgs boson search J2] can be safely used 
for this analysis as well for low mass resonances, where the bulk of the photon pj distribution 
is close to that of the photons coming from decays of the SM-like Higgs boson. This is not the 
case for high-mass resonances. Therefore a normalization uncertainty of 5% per photon pair 
is also included to account for the differences in the pj spectra of the signal photons and the 
electrons from Z —> e 1 e used to estimate the uncertainties C40| . 

The use of the BW model has an accuracy of the order \\/nix EH- A global normalization 
uncertainty equal to Tx/ is added to account for any uncertainty due to the theoretical signal 
line shape. Table |2] lists all the systematic uncertainties accounted for in the analysis. 
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7 Results and interpretation 


Table 2: Summary of uncertainties that have impact on the signal strength, applicable to events 
in all classes. 


Sources of systematic uncertainty 

Uncertainty 

Per photon 

Barrel 

Endcap 

Energy resolution, Rg > 0.94 (low //, high //) 

0 . 10 %, 0 . 20 % 

0.14%, 0.06% 

Energy resolution, Rg < 0.94 (low //, high //) 

0 . 10 %, 0.18% 

0.18%, 0 . 12 % 

Photon energy scale 

0.5% 

2 % 

Photon identification efficiency 

1 . 0 % 

2 . 6 % 

Per event 

Barrel 

Endcap 

Integrated luminosity 

2 . 6 % 

2 . 6 % 

Vertex finding efficiency 

0 . 2 % 

0 . 2 % 

Trigger efficiency 

1 . 0 % 

1 . 0 % 

R 9 class migration 

2.3% 

5.5% 

Additional normalization uncertainty 

5% 

5% 

Breit-Wigner model 

0 . 01 - 10 % 

0 . 01 - 10 % 


7 Results and interpretation 

The invariant mass spectra show no clear evidence for the presence of a new particle decay¬ 
ing to two photons. Exclusion limits are therefore computed. A modified frequentist CL S 
method 132] is employed, with the asymptotic approximation for the statistic test as described 
in Ref. Ii4l1 . Model-independent results are presented for a spin-0 and spin-2 resonance pro¬ 
duced via gluon-gluon fusion. Figure [ 6 ] shows the 95% confidence level (CL) exclusion limits 
on the production cross section times branching fraction ( aB ), obtained by combining all four 
event classes, as a function of mass for a narrow (Tx = 0.1 GeV) spin-2 resonance. Figure [7] 
shows the 95% CL combined limits for two width hypotheses, Tx = 0.1 GeV and Tx = 0.1 nix, 
as a function of mass for the spin-0 model. The results shown in Figs. [ 6 ] and [7] (left) are similar 
because of the small differences in the efficiency times the acceptance in the two different spin 
scenarios. Figure [H] shows the dependence of the limit on Tx for two values of the resonance 
mass (150 and 840 GeV) in the spin-0 model. Figure [9] shows the 95% CL exclusion limits on 
the production cross section times branching fraction as a function of the resonance mass and 
width in the spin-0 model. The expected limits lie between 6 x 10 4 pb and 4 x 10 2 pb over 
the full mass range analyzed. The observed limits are consistent with the expected sensitivity 
of the analysis in the no signal hypothesis. The largest excess is observed at ni\ ~ 580 GeV with 
a local significance of less than 2.5cr. 

7.1 Interpretation in two-Higgs-doublet model 

In this section the model-independent limits obtained for a hypothetical heavy diphoton reso¬ 
nance are interpreted in the context of the two heavy neutral Higgs bosons, H and A, predicted 
in the 2HDM llll . In this model, the production cross sections for H and A, as well as the 
branching fractions for their decays to two photons depend on two parameters, oc and /I. The 
mixing angle between H and h is given by a., while tan /3 is the ratio of the vacuum expectation 
values of the two Higgs doublets. The 2HDM cross sections are calculated by means of the 
SusHl (42l program, and branching fractions are obtained using 2HDMC |[43l . 

Exclusion regions in the tan ft versus cos(/3 — a) plane are shown only for the diphoton decay 
of the pseudoscalar Higgs boson A; no region of the phase space can be excluded for the decay 
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19.7 fb' 1 (8 TeV) 



m x [GeV] 

Figure 6: Exclusion limit at 95% CL on the cross section times branching fraction of a new, 
narrow, spin-2 resonance decaying into two photons as a function of the resonance mass hy¬ 
pothesis, combining the four classes of events. 


19.7 fb 1 (8 TeV) 



cr 1 


T 

X 


icr : 


10 “ 


10 “' 


CMS 


19.7 fb' 1 (8 TeV) 


T 


All classes combined 
r x = 10% of m x ; spin-0 


200 300 400 500 600 700 800 

m x [GeV] 


Figure 7: Exclusion limit at 95% CL on the cross section times branching fraction of a new, 
spin-0 resonance decaying into two photons as a function of the resonance mass hypothesis, 
combining the four classes of events. The results for a narrow resonance hypothesis (Tx = 
0.1 GeV) (left) and for a wide resonance hypothesis (Tx = 0.1 m x ) (right) are shown. 



of the heavy H scalar. 

Figure [10] shows the observed and expected exclusion regions for a heavy Higgs boson A of 
mass 200 and 300 GeV for the Type 12HDM. The case where H and A are degenerate in mass is 
considered. 


In Fig. 10 the region below the curve is excluded. These contour plots are similar to those 
in HIOl . The constraints obtained in this analysis on the 2HDM parameters are complementary 
to those already set by ATLAS on heavy neutral Higgs boson production using multilepton 
final states If44l and by CMS on heavy neutral Higgs bosons production using multilepton and 
diphoton final states 
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Figure 8 : Exclusion limit at 95% CL on the cross section times branching fraction of a new, 
spin -0 resonance decaying into two photons as a function of the resonance width hypothesis, 
combining the four classes of events. The results for mx = 150 (840) GeV are shown left(right). 



Figure 9: Combined observed (left) and expected (right) exclusion limits at 95% CL on the 
cross section times branching fraction of a new, spin -0 resonance decaying into two photons as 
a function of the resonance mass and width hypotheses. Contours are displayed for different 
values of <jB{X —> 77 ). 


8 Summary 

A search for resonant production of two photons is performed using 19.7 fb -1 of pp collisions 
collected at sjs = 8 TeV, in the mass range 150-850 GeV. Widths of the resonance X in the range 
0.1 GeV to O.lmx are investigated. Both spin-0 and spin-2 scenarios are considered. A fit to 
the diphoton invariant mass distribution in data is performed using a parametric model for 
the signal and a background shape obtained directly from data. No evidence for a signal is 
observed, and upper exclusion limits at 95% CL are set on the production cross section times 
branching fraction. The model-independent upper limits extend over considerably wider mass 
and width ranges than in previous searches. We further interpret these limits in the context of 
the 2HDM, presenting exclusion contours in the tan ft versus cos(/3 — a) plane. This is the first 
search for heavy diphoton resonances carried out at the LHC to be interpreted in terms of the 
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Figure 10: Observed and expected 95% CL exclusion regions for gluon-fusion production of a 
heavy Higgs boson A of mass 200 GeV (left) and 300 GeV (right) in the tan /3 versus cos(/l — a) 
plane for the Type 12HDM, assuming the H boson to be degenerate in mass with A. The regions 
below the curves are excluded. 


2HDM. 
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